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Abstract 

We take interest in a reaction-diffusion system which has been recently 
proposed in as a model for the effect of a road on propagation phenomena 
arising in epidemiology and ecology. This system consists in coupling a classical 
Fisher-KPP equation in a half-plane with a line with fast diffusion accounting 
for a straight road. The effect of the line on spreading properties of solutions 
(with compactly supported initial data) was investigated in a series of works 
starting from m- We recover these earlier results in a more general spatially 
periodic framework by exhibiting a threshold for road diffusion above which the 
propagation is driven by the road and the global speed is accelerated. We also 
discuss further applications of our approach, which will rely on the construction 
of a suitable generalized principal eigenvalue, and investigate in particular the 
spreading of solutions with exponentially decaying initial data. 

Keywords. Reaction-diffusion systems, asymptotic spreading speed, generalized 
principal eigenvalues. 


1 Introduction 

The study of the large time behavior of solutions of reaction-diffusion equations 
has been motivated by a wide range of applications, in particular in ecology and 
epidemiology. Indeed, it is well understood that propagation phenomena arise from 
the combining effect of diffusion (which accounts in an ecological context for the 
motion of individuals) and reaction (reproduction of the species), thus providing a 
mathematical model for ecological invasions or epidemics |25| . 

The most classical example is the so-called Fisher-KPP equation 

dtv = dAv + f (v), t > 0, X e M^, (IT) 

where d > 0 and the nonlinearity / S oo)), 0 < r < 1, satisfies 

/(O) = /(I) = 0, u /Al is decreasing, and —oo as v ^ -boo. (1.2) 
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The typical example is f{v) = v{l — v) and was investigated in the seminal papers |16[ 
I19| . In particular, this equation is well-known to admit positive traveling waves, 
namely planar solutions moving with a constant shape and speed through the domain 
and connecting the two constant steady states 0 and 1 dKEKIl]. The minimal wave 
speed, which can be explicitly computed as 

^*KPP — 2\/ df'{0), 

is also the asymptotic spreading speed of solutions of the Cauchy problem (HU 
with compactly supported initial data [2]. The striking feature of these results is 
that, under the KPP assumption (|1.2I) . the speed ^KPP is fully determined by the 
linearization of (|l.ip around u = 0. 

Various extensions of the KPP equation have been investigated in the past years. 
In particular, a tremendous effort has been dedicated to the inclusion of heterogene¬ 
ity, which is a common yet mathematically challenging feature of more realistic eco¬ 
logical models |25) . Among the vast literature, we refer in particular to [a El [710 [27] 
for generalizations of the results mentioned above in the spatially periodic framework, 
in which we will also place ourselves below. Let us also mention that further hetero¬ 
geneous problems, as well as other types of nonlinearities have also been studied: we 
refer again to 0 and the references therein, as well as to the thorough reviews 

More recently, Berestycki, Roquejoffre and Rossi proposed a new model specif¬ 
ically devised for studying the role of roads in biological invasions m This issue 
was motivated by the empirical observation that not only diseases, but also various 
ecological species, tend to spread faster along roads and other transport lines. 

Their model is a system of two parabolic equations, whose unknowns u{t,x) and 
v{t,x,y) should be interpreted as the density of a same population but on different 
domains, respectively the road and the field: 

{ dtu — Dd^u = i'{x)v{t, X, 0) — fi{x)u, t > 0, x S M, 

dtv - dAv = f{v), t > 0, (x, ?/) G n, (1.3) 

—ddyv{t, X, 0) = fi{x)u — iy{x)v{t, x, 0), t > 0, x G M. 

From now on and consistently with this interpretation, we will refer to the upper 
half plane 12 = {(x, y) G : y > 0} as the field, and to the line {(x, 0) : x G M} 
as the road. 

In |llj as well as throughout this work, the field equation (second line of (ll.3p f is 
still assumed to be of the KPP type in the sense that the function / satisfies (II.2p . 
However, it is coupled with the road equation (first line of (|1.3I) 1 through a Robin 
type boundary condition at y = 0 (last line of (II.3h ). Roughly speaking, fiu is the 
proportion of individuals jumping from the road to the field, and uv the proportion 
of individuals jumping from the field to the road. Note that the field boundary 
condition is normalized so that, when no reproduction takes place (/ = 0), then the 
total mass of the population is conserved (we refer to Ha for a complete argument). 

When /i and p are positive constants, Berestycki, Roquejoffre and Rossi Ha 
proved that the solution (with compactly supported initial data) spreads along the 
road with some speed c*{D) > c’^pp = df'{0) and that, moreover, this inequality 

is strict if and only if D > 2d. This in particular means that, to accelerate the propa¬ 
gation, fast diffusion on the road is enough even though individuals do not reproduce 
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there. One may also note that this diffusion threshold D = 2d for acceleration does 
not depend on the values of the constants /i and u, an observation which will be 
strengthened by our results. 

A more general formula for the diffusion threshold, including the case when 
transport or reproduction take place on the road, was given in Uni: in particular, 
when reproduction rates are identical on the road and in the field, then acceleration 
occurs as soon as D > d. Related results were also obtained in various frameworks 
such as fractional diffusion [4], ignition type nonlinearities [H [E] and nonlocal 
exchange |23| . The issue of the spreading speed in oblique directions away from the 
road was investigated in mm and revealed interesting properties where the shape 
of the propagation is no longer dictated by a Huygens principle, and may not even 
be convex. Lastly, we refer to |26| where the restriction of the KPP problem (II.3p 
to a truncated field, with a Dirichlet boundary condition and constant exchange 
coefficients, was studied. It turns out that such a truncation plays an essential role 
in our arguments. 

While all the above papers dealt with systems that are invariant in the x-direction 
(let us point out here that the presence of the road itself entails strong heterogeneity 
with respect to the y-direction), our work will stand in a more general spatially 
periodic framework where the results of |11| (as well as some of [101 I26| ) will be 
extended. Even though the heterogeneity will be limited to the exchange terms, we 
believe our arguments could be extended to more general periodic systems and will 
briefly discuss this in the last section. 

1.1 Main results 

From now on, we consider the system (11.311 where D, d > 0, the function / satisfies 
the KPP assumption (|1.2I) . and /r(x),z^(x) are L-periodic exchange coefficients in 
C^’'’(M) such that 

/r 0, u 0. 

As mentioned above the main difference with the original model of m lies in the 
x-dependence of the exchange terms, which is motivated by the ubiquity of het¬ 
erogeneity in the applications. The periodicity provides a suitable mathematical 
framework, while still requiring a new approach avoiding explicit computations. We 
also point out that /r and u may occasionally be 0, which of course could not hap¬ 
pen in the homogeneous exchange case. This can be interpreted as the presence of 
"walls" blocking locally the motion of individuals between the field and the road 
(but not necessarily both ways simultaneously). 

For convenience, we introduce the notations 

fiQ := min/i < /i(x) < /ri := max/i, i'q := mini/ < i/(x) < vi := maxi/. 

The system (11.31) is supplemented with bounded and non-negative initial data 

r u\t=o = uo in M, , . 

\ u|t=o = xo in 11. ^ 

As will be shown in Section [2l this problem is well-posed under generic conditions 
on the initial data. 
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In this paper, we take interest in the large time behavior of the solution and 
more specifically in its spreading speed in the x-direction or, in other words, along 
the road. By spreading, we mean here that the solution should converge to a positive 
steady state on some expanding in time domains. Therefore, we will begin with a 
Liouville type result to describe such positive steady states: 

Theorem 1.1. There exists a unique positive and bounded stationary solution {U, V). 
This solution is L-periodic in x, has a positive infimum and satisfies y(x,+oo) = 1 
where the limit is uniform with respect to x G M. 

Unlike in the homogeneous exchange case, it is not possible to explicitly compute 
this positive steady state and, therefore, both its existence and uniqueness will be 
inferred from the (linear) instability of 0 and the KPP assumption. The proof, which 
will be performed in Section 12.21 relies on the fact that the road-field system satisfies 
a comparison principle, namely Propositions 12.11 and 12.21 Indeed, although it is 
non-standard due to the 1D-2D coupling through a boundary condition, this model 
shares many features with monotone systems. We refer the reader to |21| and the 
references therein for a study of abstract monotone systems. 

From the proof of Theorem II. H it will be clear that this steady state is stable. One 
may thus expect it to attract solutions of (II. leading to the spreading dynamics 
that we wish to investigate: 

Theorem 1.2. There exists c*{D) > c*j^pp = 2Y'df'(()) > 0 such that any solution 
{u,v) of ([I3])-([I31) with non-negative, continuous and compactly supported initial 
data {uo,Vo) ^ (0,0) satisfies: 

• for all c > c*{D) and R> 0, 

lim sup {u{t,x) + v{t,x,y)) = 0] 

t->-+oo x<-ct , 0<y<R 

• for all 0 < c < c*{D) and R > 0, 

lim sup {\u{t, x) — U{x)\ + \v{t, x,y) — V{x,y)\) = 0. 

t^-Coo o>x>-ct , 0<y<R 

The above statement deals with spreading in the left direction. As the symmetric 
problem x to- — x satisfies the same assumptions, there clearly also exists a spreading 
speed in the right direction. However, the left and right spreading speeds may be 
different. 

The proof is directly inspired by the single (spatially periodic) equation. It relies 
on the construction of a family of principal eigenvalues denoted by A(q:), which arises 
when looking for exponential solutions 

e"(^+^^)(U,U) 

of the linearized problem around the invaded unstable state (« = 0, u = 0), where 
a > 0 and U{x), V{x,y) are positive and L-periodic with respect to x. Due to the 
unboundedness of the domain, the usual notion of principal eigenvalue is not uniquely 
defined and there exist several definitions of "generalized" principal eigenvalues |13| . 
Our generalized principal eigenvalue A (a) (see Theorem 13.11 for its main properties) 
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will be constructed in Section [3] by truncating the field in the y-direction and passing 
to the limit to the whole road-field domain. 

Then in Section HI we will prove Theorem 11.21 and show that the spreading speed 
c*{D) is characterized by the following formula: 


c*{D) 


min 

q ;>0 


—A(a) 
a 


This is the natural extension of a similar formula for the minimal wave speed of the 
single spatially periodic equation [3 El, which in the homogeneous case (HID easily 
reduces to 


Cft'pp = mm ■ 
a>0 


da^ + f'{0) 


a 


(1.5) 


Therefore, whether the road accelerates the propagation or not depends on how 
—A(a;) compares to + f'{0) (note that it is already stated in Theorem 11.21 that 
c*{D) > c*j^pp)- In fact we will prove that — A(q:) is strictly larger than da'^ + f'{^) 

if and only ii D > d and a > see Proposition 15.11 Recalling that a denotes 

the exponential decay of the ansatz as x —)• — oo, the formal argument reads as 
follows: when looking at the linearized problem (13.3D later on, the intrinsic growth 
rate of the road population is Da^ (due to motion of road individuals in the original 
problem) and, in the field, it is da^ f'{0) (due to both the motion of individuals in 
the field and their reproduction). Thus, one may indeed expect the road to lead the 


propagation if and only \i a > ^J Notice then that, as the exponential decay 

of the KPP traveling wave with minimal speed is \J ^ which is also where the 
minimum in (II. 5D is reached, the diffusion threshold D = 2d immediately arises. 

The above argument will be made rigorous in Section [5l leading to the following 
theorem: 


Theorem 1.3. The strict inequality c*{D) > ^KPP 
Furthermore, 

c*{D) 


holds if and only if D > 2d. 


0 < lim inf ■ 
D^oo 


yfD 


, c*{D) 

< hm sup — 1 =^ 

D^OO yD 


< - 1 - 00 . 


We highlight the fact that we recover the exact same diffusion threshold D = 2d 
as in m, not only in spite of the heterogeneous exchange, but also in spite of the 
fact that fj, and n may occasionally be 0. 


The above discussion also leads one to expect that if the initial data are slowly 
exponentially decaying as x —>■ —oo, then the solution may no longer be accelerated 
by the road. This can be easily interpreted since the effect of diffusion gets weaker 
as the solution is relatively flat. More generally we will prove in Section [5.31 that, as 
in the single equation case, the family of eigenvalues —A(a) also characterizes the 
spreading speed of solutions of (|1.3p with any exponential decay: 

Proposition 1.4. Let a* > 0 be the smallest solution ofc*{D)a = —A{a), and {u, v) 
be the solution of (fr3D-(fL4l) with non-negative, bounded and continuous initial data 
{uo,Vo) satisfying also 

( m(0)e"* < uo{x) < Me“^, 

^ ' I m{y)e°^ < vo{x,y) < Me“^, 
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for some 0<a<a*,M>0, and m a continuous and positive function of y ^ 
[0,oo). 

Then {u,v) spreads along the road (in the same sense as in Theorem \1.2\) with 
speed 


In particular, if either D < d or D > d and a < then c{a) = da + is 

also the speed of solutions of dm with exponential decay of order e“^. 


Combining Theorem 11.21 and Proposition 11.41 it is clear by a simple comparison 
argument that any solution which decays at least as as x ^ — oo at time t = 0 
spreads with the speed c*{D). In particular, it turns out that the behavior of the 
initial data in the y-direction barely matters, which is in some sense natural but also 
not trivial as the eigenfunctions associated with —A(a) may decay as y —)• +oo. 


Outline of the paper. The paper is organized as follows: in Section [2] we inves¬ 
tigate well-posedness of the Cauchy problem (II.3p and its stationary solutions. In 
Section [3] we construct the principal eigenvalue of the linearized problem and the 
associated exponential solutions. Section |4] contains the proof of the main spreading 
Theorem 11.21 We investigate in Section [5] the acceleration properties and the limit as 
H —>• oo of large diffusion on the road, that is Theorem 11.31 as well as the spreading 
properties for exponentially decaying initial data. 

Lastly we will discuss in Section [6] some possible extensions of our work. We will 
focus in particular on the truncated field case, i.e. the restriction of (|1.3I) to the 
domain {0 < y < R}, and on the case when reproduction also takes place on the 
road. We will see that most of our results can be extended in both those frameworks 
with only minor modifications. We will then discuss the general spatially periodic 
problem and point out some remaining open questions. 


2 The Cauchy problem and stationary solutions 

2.1 Preliminaries: the evolution problem 

We begin by noting that the parabolic system satisfies the comparison principle. 
From now on, we say that {u, v) is a subsolution (respectively supersolution) of (11.31) 
if it is satisfied with inequalities < (respectively >) instead of equalities, and if both 
functions u and v are continuous. 

Proposition 2.1. Let {u,v) and {u,v) be respectively a subsolution bounded from 
above and a supersolution bounded from below such that {u,v) < {u,v) at time t = 0. 

Then, either {u,v) < {u,v) for all t > 0, or there exists T > 0 such that {u,v) = 
{u,v) for all t <T. 

Proof. The proof is the same as in m and we omit the details. The argument relies 
on a combination of the strong parabolic maximum principle and parabolic Hopf 
lemma, together with the key assumption that y and u are non-negative. □ 

This comparison principle also extends to generalized sub and supersolutions. 
We again omit the proof and refer to Proposition 3.3 in m for the details. We write 
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their result below for the convenience of the reader and because we will use it several 
times throughout this paper: 

Proposition 2.2. Let E C (0, +oo) x M and F C (0, +oo) xLl be two open sets. Let 
{ui,vi) and (^ 2 ,^ 2 ) be two subsolutions of (II.3p . bounded from above and satisfying 

ui < U 2 on dE n ((0, + 00 ) X M), vi < V 2 on dF n ((0, + 00 ) x Ll), 


in E, 
otherwise, 

in F, 
otherwise. 

If it satisfies 


and define the functions 

u{t, x) : = 

v{t,x) := 


max{ui, U2} 

U2 

max{ui,U2} 

V2 


yft, x) > U 2 {t, x) ^ v{t, X, 0) > vi{t, X, 0), 
v{t, X, 0) > V 2 {t, X, 0) u{t, x) > ui{t, x), 

then, any supersolution {u,v) of dLl such that u<u and v_<v at time t = 0, also 
satisfies u<u and v<v for all time t > 0. 

Thanks to those comparison principles, we can prove as announced that the 
Cauchy problem is well-posed. 

Theorem 2.3. The Cauchy problem (ll.3p - (ll.4p . where the initial datum is bounded, 
non-negative and Holder continuous, has a unique global non-negative and bounded 
solution. 


Proof. Note that uniqueness immediately follows from the previous proposition, so 
it only remains to prove the existence part. The proof is again very similar to m, 
by constructing in the spirit of |24) a monotonic sequence of solutions (u”,u”)n>i to 


dtu"'— Dd'^u^ + n{x)u"'= u{x)v^ ^{t,x,0), t>0,x£R, 

u^\t=o = uo, x £ M, 


( 2 . 1 ) 


and 


dtv^ — dAv^ = /(u"'), t > 0, {x, y) G Ll, 

< —ddyv"'{t,x,0) = p,{x)u'^ — i'{x)v'^{t,x,0), t>0,x£M., (2-2) 

, v"‘\t=o = vo, X £ M, 

where > 0 will be suitably initialized when n = 1 in (|2.1I) later on. 

By classical parabolic theory |20) both problems shall be solvable at each step, 
and the minimum principle ensures that u”, u” stay non-negative. In order to pass to 
the limit and obtain a solution of (II.3p . we need some L°° estimates on (u”, u”) uni¬ 
formly with respect to n. This is done by exhibiting a bounded supersolution (which 
will also serve in the initialization of the sequence) and applying the comparison 
principle. 
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First let Uq{x) be the (positive) principal eigenfunction of the periodic problem 

f -DUq + fi{x)Uo = XqUq, X e M, 

1 Uq is L-periodic. 

Note that from our assumptions on p, it is clear that this principal eigenvalue Aq > 0. 

Now let 

{U{x),V{x,y)) :=K(Uoix),C{l + e-^y)) (2.3) 

for some parameters K,C,oj > 0 to be adjusted shortly. Fix first C < and 

then (jj > gy KPP hypothesis (jl.2p we can next choose K > 0 large 

enough such that < —duj^, and some straightforward computations show that 

(U,V) satisfies 

( -DU" + fi{x)U >iy{x)V{x,0), 

\ -dAV>f{V), 

[ —ddyV{x,0) + u{x)V{x,0)>fi{x)U. 

Note that this positive stationary supersolution is bounded from above and away 
from zero: up to increasing K > 0 (for hxed C, oj) we can therefore assume that the 
initial data 0 < (uo,vo) < {U, V). 

Initializing v^{t,x) = V(x) in (12.11) . a straightforward inductive application of 
the comparison principle, Proposition 12.11 shows that 

U> ...>u"> > ... > 0 and V > ... > v" > > ... > 0 

for all n > 1, and therefore 

Vt > 0, X € M, y > 0 : x, y) \ u{t, x, y) and v"{t, x) \ v{t, x) 

as n —)• oo pointwise. 

In order to show that the pair {u, v) satisfies (jl.3j) together with the initial con¬ 
dition m, we claim now that the sequence {u”,u"'}n>i is relatively compact in 
the local uniform, the C'|^^((0,oo) x 11), and the Cj^‘^{{0,oo) x H) topologies. More 
precisely, for any i? > 0 and T > 0, the uniform boundedness of implies uniform 
VFp’^((0,T) X Bji) estimates on the sequence u" for any large p, and in particular 
it is locally uniformly Holder continuous by Morrey’s inequality. Then by parabolic 
estimates m on the field equation the sequence u” is locally uniformly Holder con¬ 
tinuous, as well as its spatial derivatives of order 1 for all t > 0 and up to the 
road boundary {y = 0}. Furthermore, by standard interior Schauder estimates, its 
time derivative and space derivatives up to the order 2 are locally uniformly Holder 
continuous inside the field. This is enough to obtained the desired compactness. 

Clearly any cluster point must agree with the previous pointwise limit {u",v'") 
{u,v), so by standard uniqueness and separation arguments we conclude that the 
whole sequence converges in the strong topologies. This is enough to pass to the limit 
in and retrieve (|1.3p for all positive times. Moreover, (n, u) = 

is continuous up to time t = 0, and the proof is achieved. □ 





2.2 The Liouville type result 

We look now for stationary solutions, i.e. solutions {U, V) of 

—DAU = {x, 0) — fj,{x)U, a; G M, 

< -dAV = f{V), {x,y)en, (2.4) 

—ddyV{x,0)=fj,{x)U{x) — u{x)V{x,0), x G M, 

and prove Theorem 11.11 Observe that, in the statement of Theorem m the peri¬ 
odicity of {U,V) immediately follows from uniqueness since any (lattice) translation 
of {U, V) is again a solution. Moreover, the limit V{x, -|-oo) = 1 is a straightforward 
consequence of the positivity of the infimum: indeed, up to extraction of a subse¬ 
quence, V converges as y —)■ -|-oo to a positive and bounded stationary solution of 
—dAV = f{V), which may only be 1 by the KPP assumption. 


Let us now begin the proof of Theorem 1 1.1 1 with the existence part: 


Lemma 2.4. There exists at least one non-trivial, non-negative and bounded solution 
of (ITD. 

Proof. We construct below some non-negative subsolution {U_,V_) of (|1.3p which is 
bounded and non-trivial. We use the classical method and hrst let ii > 0 large 
enough such that the principal eigenvalue X{R) of —dA in the ball i?/j(0,0) C 
with Dirichlet boundary condition is less than /'(0)/2. If $jj(x,y) > 0 is the corre¬ 
sponding principal eigenfunction, set 


U{x):=0, V{x,y) 


^Rix,y-R-l) ii {x,y) £ Br{0,R-^1), 
0 elsewhere. 


Choosing then a constant C > 0 small enough and thanks to the regularity of /, 
{U_,V_) = C{U,V) will automatically be a subsolution of (12.4p . 

By Theorem 12.31 there exists a solution {u,v) of ()1.3p with initial datum 
Furthermore, as (II.3p satisfies a comparison principle (see Propositions 12.1 1 and 12.21) . 
it lies above {U_,y_) for all time and, furthermore, it is non-decreasing with re¬ 
spect to time. On the other hand, up to decreasing C, we can also assume that 
{U_,V_) < {U,V) where {U,V) is the supersolution (j2.3p constructed in the pre¬ 
vious section. In particular, {u, v) is bounded from above, hence it converges as 
t —>• -|-oo. Recalling that this supersolution is uniformly bounded, one can use stan¬ 
dard parabolic estimates (see also the end of the proof of Theorem 12.3p to get that 
{u,v) converges locally uniformly to a stationary solution {U,V) of (12.4p . It is clear 
by construction that {U, V) is non-trivial, non-negative and bounded. □ 


Lemma 2.5. Let {U,V) be a non-negative, non-trivial, and bounded solution of 
UM- then infft [/ > 0 and inf^P > 0. 

Proof. We show first that P > 0 is bounded away from zero, and positivity of U will 
follow. 

If P >0, then the classical strong maximum principle, applied in the field, 
implies P > 0. Using again the first Dirichlet eigenfunction of —dA in a large ball 
Br C Ll and moving the ball it is easy to show as in m Proposition 4.1] that 


Vr > 0, 


(2.5) 
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inf P > 0. 

y>r 



In other words, V is uniformly positive far from the road. Arguing by contradiction, 
assume that V (x„, y„) \ 0 for some sequence (x„, yn)n>o such that \ 0. Shifting 

Unix) :=Uix + Xn), Vnix, y) ■.= Vix + Xn, V + yn), 

Unix) ■■= yix + Xn), l^nix) := l^ix + Xn), 

elliptic estimates m and periodicity of y and show, up to extraction of a sub¬ 
sequence, that Un, Vn, Hn, and Vn converge locally uniformly to some U,V,jl,D 
satisfying the same system in the same domains, and 

V >0 in n, F(0,0) = 0. 

In fact (12.5p implies that V > 0 in 12. Since —dAV = fiV) > 0 in some neighborhood 
of (0,0), Hopf lemma finally yields the following contradiction 

0 > -ddyViO, 0) = yi0)Ui0) - i>(0)y(0,0) > 0. 

Therefore, inf-^V > 0 as announced. 

Assume now there exists ixn)n>o such that U{xn) \( 0. As above, up to trans¬ 
lating in X and up to extraction of a subsequence, we may assume that 

Uix) >0 in M, 17(0) = 0. 

Here, abusing of notations, we still denote by U, V, y and u the shifted functions 
satisfying the same system. Since V is positive, we have 

-DAU + yU = uV{-,Q) > 0. 

Hence, by the strong maximum principle, we get that 17 = 0. But then, i'i-)V (-,0) = 
0, and choosing any xq such that z^(xo) > 0 contradicts the fact that V has positive 
infimum. This completes the proof. □ 

As mentioned above, it follows from the positive inhmum of V and the KPP 
assumption that 

V{-, -l-oo) = 1. 

As explained at the beginning of this section, it only remains to prove the uniqueness 
of positive and bounded stationary solutions. 

Lemma 2.6. There exists at most one non-negative and non-trivial bounded solution 
to (IO) . 

Proof. The idea is similar to |8] (we refer more specifically to the proof of their The¬ 
orem 2.4). The argument strongly relies on the uniform positivity from Lemma [2.51 
which followed from the instability of (0,0) as a steady state of (11.31) . As this insta¬ 
bility may be reinterpreted in terms of the generalized principal eigenvalue defined 
in Section [3l we also refer to the deeply related Liouville type results of |9] for the 
general heterogeneous KPP equation. 

Assume (t7i,Vi) and {U 2 ,V 2 ) are two such solutions, and let 

r:=sup{0>O : iUi,Vi)>9iU2,V2)}. 
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Since both solutions are bounded away from zero (Lemma 12.51) and from above, 6* is 
finite and positive. We show below that 6* >\ and therefore (C/i, Vi) > {U 2 , V 2 ). By 
symmetry we would also retrieve the opposite inequality, hence [Ui,Vi) = {U 2 ,V 2 ). 
We proceed by contradiction and assume that 6* <1. By continuity there holds 

P{x) := Ui{x) - e*U 2 {x) > 0, Q{x,y) := Vi{x,y) - B*V 2 {x,y) > 0, 

and by definition of 6* there exists a sequence {xn,yn)n>o such that P{xn) 0 or 
Qixn,yn) \0 (or both). 

We first consider the second case. Note that 


Q{x,y) 1 - 0* > 0 


as y ^ + 00 , thus we may assume without loss of generality that the sequence y^ is 
bounded and converges to some ?/oo > 0. 

Assuming first that ?/oo > 0, shifting in the x direction as before, and extracting 
a subsequence we can assume that Q{0, yoa) = 0 and Q >0 satisfies 

-dAQ = f{Vi) - e*fiV2) > fiVi) - f{e*V2). 


Here we took advantage of the KPP hypothesis {f{v)/v decreasing) and of the in¬ 
equality 6* < 1. Thus 

- dAQ + aQ>0 (2.6) 


in M X (0, + 00 ), where 


a{x,y) = 


fjVi) - f{0*V2) 

Vi - e*V2 


is bounded uniformly in x,y because Vi, V 2 are and / is Lipschitz. Since Q attains an 
interior minimum point (0,yoo) S the strong maximum principle implies Q = 0, 
which contradicts the strict inequality in (|2.6I) . 

Assuming now that y^o = 0 and shifting again, we may take without loss of 
generality <5(0,0) = 0 and, according to the previous argument, <5 > 0 in n. Since 
(12.6p still holds, Hopf lemma shows that dyQ{0,0) > 0, which is impossible since 


0 > -ddyQ{0,0) = ^^{0)P{0) - i/(0)Q(0,0) = yi{0)P{0) 


and P > 0 by definition of 9*. 

Since we just proved that <5 > 0 up to the boundary we only have to prove 
that P > 0, which will yield the desired contradiction. If not, then shifting again, 
extracting a subsequence, and with the same abuse of notations, we may assume 
that P(0) = 0, where P still satisfies the first equation 


—DAP + fi{x)P = i'{x)Q{x,0), X G M. 

Since <5 > 0 and y,(x),u{x) > 0 the strong maximum principle implies P{x) = 
P(0) = 0. The previous equality immediately implies 0) = 0, thus Q{xq, 0) = 

0 at any point xq such that ^{xq) > 0. According to the previous argument this is 
a contradiction, and the proof is complete. □ 
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3 A generalized principal eigenvalue 


In this section, we take interest in the linearized problem of ()1.3p : 


dtu — DAu = X, 0) — fi{x)u, 

dtv — dAv = f'{0)v, 


t > 0, X € M, 
t > 0, (x, y) e n, 


(3.1) 


—ddyv{t,x,0)=y{x)u — iy{x)v{t,x,0), t>0,x£ 


We will seek for exponential solutions of (13.11) moving with constant speed, i.e. of 
the form 

{u, v) = {U{x), V{x, y)), c> 0, a > 0 (3.2) 

for some x-periodic and positive U, V. By analogy with the single equation, we ex¬ 
pect that there exists a critical speed c*{D) such that these exponential pulsating 
wave solutions exist if and only if c > c*{D). Furthermore, we expect this critical 
speed c*{D) to be also the spreading speed of solutions of (|1.3p with compactly sup¬ 
ported initial data. 

Fixing a parameter a > 0, we will denote for convenience 

Li{U, V) := -DU" - 2DaU' + {-Da^ + y{x))U - iy(x)V(x, 0), 

L2iU, V) := -dAV - 2dad^V + {-da‘^ - f{0))V, 

E{U, V) := -ddyV(x, 0) + v{x)V{x, 0) - y{x)U{x). 

Then, plugging the previous ansatz into (|3.ip leads to the following eigenvalue prob¬ 
lem 

Li{Ua,Va) = AUa in M, 

L2{U^,Vo,) = AV^ 


. £;(!/«, W.) = 0 

together with the conditions 


m 

in 


X (0, -l-oo), 


(3.3) 


and 


Ua is L-periodic and positive, 
Va is L-periodic and positive, 

A = —ac. 


In this section, we will deal with the well-posedness of this eigenvalue problem. The 
main difficulty comes from the unboundedness of the domain, which typically entails 
the non-uniqueness of the eigenvalue associated with positive eigenfunctions; thus 
the need of some well-chosen "generalized" eigenvalues. We refer to |13j and the 
references therein for many results on such eigenvalues and their applications (for 
single equations). 

Here, we will construct a generalized principal eigenvalue by first truncating the 
domain in the y-direction, then passing to the limit. This will lead us to the following 
theorem. 

Theorem 3.1. There exists a concave (hence continuous) function a > 0 i— )• A(a) 
such that : 
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(i) the system with A = A{a) admits a positive and L-periodic solution 

{Ua,Va) such that, up to some normalization, 

Uaix)<l, Vaix,y) <C{l + y), 

for some C > 0 (possibly depending on a); 

(ii) the system (13.31) admits no positive and L-periodic solution for any A > A(a). 

The rest of this section is devoted to the proof of Theorem 13.11 

3.1 Principal eigenvalue on truncated fields 

Here and henceforth, a > 0 is a fixed parameter. As announced, we first seek for a 
principal eigenvalue in the truncated fields 

Lin := T X (0, R), 

where T denotes with some slight abuse of notations the one dimensional torus 
More precisely, we seek Aij(Q:) G M such that there exists 

{Ua,R, Va,R) > 0 periodic with respect to x 

solution of the following system: 

Li{Ua,R, Va,R) = AR{a)Ua,R in T, 

< L 2 {Ua,R,Va,R) = AR{a)Va,R in P/J, (3.4) 

^ E{Ua,R, Va,R) = 0 = Va,R\y=R in T. 

The compactness of the domain, together with the cooperative nature of the road- 
field system, now provides a much more convenient mathematical framework where 
the classical Krein-Rutman theory applies. However, because of the non-standard 
coupling through a boundary condition, and anticipating the need to pass to the 
limit as ii —)• -|-oo, we include the details of the proof. 

Proposition 3.2. There exists Mq, such that, for any M > M„ and given any 
functions gi G and §2 G there exists a unique solution U,V to 

' Li{U,V) + MU = gi inT, 

< L2{U,V) + MV = g2 inLlR, (3.5) 

E{U,V) = t) = V\y=R inT, 

with U G (72’'-(T) and V G C^’’’(nR) n C^’’’(Il^). 

Moreover, (1^ enjoys a strong maximum principle: gi,g 2 > 0 implies U,V >0, 
and if either gi ^ 0 or g 2 ^ 0 then U > 0 inT and V >0 m T X [0, i?). 

Proof. Our proof is divided into three steps. The first step consists in constructing 
families of sub and supersolutions, which will serve in proving both the maximum 
principle (step 2), and the existence of a solution (step 3). 

Step 1: construction of barriers. We first seek positive supersolutions of the form 

{U{x),V{x,y))=K {1,0(1+ e-‘^y)) 
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for some large K > Q. Choosing for example w = 1, C = /ri/d, an explicit compu¬ 
tation shows that if 

M > M„ := max ^d{a^ + 1) + /'(O), Da^ + 2^^} (3.6) 

then we can choose K > 0 large enough (depending only on maxg^i and max <72 if 
positive) such that the pair (?7,1/) > 0 is indeed a supersolution of (13.51) . By linearity 
and up to increasing K (depending on min gi and min g 2 if negative), ([/, = 

-(u,v )<0 is also a negative subsolution. Moreover, it is easy to check that these 
are strict sub and supersolutions, i.e. they are not solutions. 

Step 2: strong maximum principle and uniqueness. For fixed data gi,g 2 > 0, let 
{U,V) be any solution of (13.5p . One can then easily check that {U_,y_), as defined 
above, is a strict subsolution of (13.5p for any K > 0. Now assume by contradiction 
that U or V take negative values somewhere in their respective domain. Because 
the previous subsolutions (U^, V_) are negative up to the boundary we can define 

e* = min{0 > 0 : {U,V) > 9{U,V)} G (0,+oo). 

Then (U — V — 9*V_) > 0 is a supersolution of (13.5p and attains zero either on 
the road or in the field. In the former case, the strong maximum principle implies 
that U = 0*U_, which in turn implies that V{y = 0) = 6*V_{y = 0) and, thanks to 
the Hopf lemma, V = 9*V_. This contradicts the fact that {U, F) is not solution of 
(13.5p . A similar argument leads to the same contradiction in the latter case, and we 
conclude that ([/, V) > 0. 

The same argument as above shows that {U,V) > 0 as soon as 51 ^ 0 or 52 ^ 0 
and we omit the details. Moreover, by the linearity of (13.5p . it immediately follows 
that ( 0 , 0 ) is the unique solution when gi = 0 and 5^2 = 0 and, therefore, the 
system (13.51) always admits at most one solution. 

Step 3: existence. Arguing as in the proof of Theorem 12.31 it is easy to see that 
the linear evolution problem 

dtu =—{Li + M){u,v) + gi in T x (0,-|-oo), 
dtv =-{L 2 +M){u,v) + g 2 in X (0,-|-oo), 

E{u, v) = 0 = v\y=ji in T X (0, -|-oo), 

, {u,v)t=o = {U,V), 

is solvable. Since the initial data {U_,y_) = —{U, V) is a subsolution of the stationary 
problem the solution (u, v) is non-decreasing in time, and by the maximum principle 

Vt>0: {LL,V) <{u,v) <{U,V). 

Letting t —)• 00 we get by standard parabolic estimates and monotonicity that (rt, v) 
converges to a stationary solution {U,V) of (13.5p . Standard elliptic estimates give 
the desired regularity for U{x),V{x,y) and the proof is complete. □ 

According to Proposition 13.21 above, the system (13.5p is well-posed for M large 
enough. By the Krein-Rutman theory, we will get as a consequence that: 
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Proposition 3.3. For any R > 0 and a > 0 there exists a unique principal eigen¬ 
value Kii{a) of (j3.4jl associated with positive and periodic eigenfunctions {Ua,R,Va,R) 
of ()3.4I) . which are also unique up to multiplication by a positive factor. 

Furthermore, — A/j(q;) < Ma, where Ma is defined by Proposition \d.S\ and more 
precisely in dM]). 

Proof. Let us fix any 0 < r' < r. Define first the Banach space 

X := {{U,V) G C^’"'(T) X x [0,R]) : V\y=R = 0} 

with its natural Holder norm, and the positive cone 

:= {{U,V) G C7i’^'(T) x x [0,R]) : U,V>0, V\y=R = 0}. 

Note that has non-empty interior in X (take for instance U (x) = 1 and V(x, y) = 
l-y/R). 

Choose now any M > Ma'. by Proposition 13.21 one can define the operator 
T '. X X, which to any given pair ( 51 , 52 ) associates the unique solution (U,V) 
of (|3.5p . Note that T is compact. This follows from elliptic estimates m and the 
fact that, for any functions 51 and 52 , the L°° norm of (U, V) = T( 5 i, 52 ) is bounded 
from above by some constant which only depends on || 5 i||oo and II 52 II 00 (see the proof 
of Proposition 13.21 above). Moreover, using again Proposition 13.21 T{X~^) C X'^, 
and T is even strongly positive in the sense that : if ( 51 , 52 ) G X^ \ {(0,0)} and 
{U,V) = r( 5 i, 52 ), then {U,V) > 0 for 5 G [0,42) as well as dyV{x,R) < 0 (by Hopf 
lemma), which implies that {U,V) is an interior point of X^. 

By the Krein-Rutman theorem we conclude that there exists a positive eigenvalue 
(Tm > 0 of T associated with positive (and bounded) eigenfunctions Um,Vm > 0. 
Transferring as usual in terms of the original PDE we have therefore an eigenvalue 
^ii(ct) = ~ of (13.4h . associated with the same positive and periodic eigen¬ 

functions. 

It now remains to check that such a Ar^o) is unique (in particular, it did not 
depend on the choice of M > Ma above). Let Ai and A 2 be two eigenvalues of 
(j3.4p associated with positive and periodic eigenfunctions, respectively {Ui,Vi) and 
{U 2 , V 2 ), and set 


9* := min{6( > 0 : 9{Ui,Vi) > {U 2 , V 2 )}. 

li9*Vi = V 2 , then trivially Ai = A 2 (the field equation does not depend on 17). Oth¬ 
erwise, we know that 6*{Ui, Vi) > {U 2 , V 2 ) and, by Hopf lemma and x-periodicity, 

dy{e*Vl-V2)\y=R<0. 

Hence, by construction, either 9*Ui = U 2 for some xn G T, or 9*Vi = V 2 for some 
(xo,yo) eTx [0,42). 

In the former case, then 

0 > rLi([/i, Hi)(xo) - Li([/ 2 , V 2 ){xo) = (Ai - A 2 )U 2 {xo), 

hence Ai < A 2 . One gets the same conclusion in the latter case if the contact 
point occurs in the interior, namely 50 > 0. Otherwise, 0*Vi > 9*V2 in D/j and 
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6*Vi{xo, 0) = V 2 (xo, 0), which leads to a contradiction by Hopf lemma. We conclude 
that Ai < A 2 , and thus by symmetry Ai = A 2 . Using the usual combination of the 
strong maximum principle and Hopf lemma, it is now straightforward to check that 
6*{Ui,Vi) = {U 2 ,V 2 ), namely that the principal eigenfunction pair is unique up to 
multiplication by a positive factor. 

Finally, we end the proof by noting that the upper bound —A/j(a) < follows 
from the fact that AR{a) = ^ — M > — M for any M > M^. □ 

3.2 Passage to the whole field 

Before we pass to the limit as i? —)• + 00 , we need some monotonicity and bounds on 
the function R 1 —?■ A/j(a). This is the purpose of the proposition below: 

Proposition 3.4. For any a > 0, the function i? > 0 1 —)• —AR{a) is increasing. 
Moreover, 


max + f\0) — , Da^ — Aq| < —Aj^(q;) < M^, 

where was defined in Proposition \3.‘A and Aq, £ [min fi, max p] is the principal 
eigenvalue of —D-^ — 2aD-^ + l^ix) on the torus. 

Proof. Let us briefly sketch the proof of the monotonicity with respect to R. Fix 
Ri > R 2 > 0 and let Ai = Ajj^(a), A 2 = A/j 2 (q^) be associated with the positive 
and periodic eigenfunctions (17i,Vi) and {U 2 ,V 2 ). Since Ri > R 2 , one can proceed 
as before and, multiplying {Ui, Ui) by a well-chosen constant 9 > 0, assume without 
loss of generality that {Ui,Vi) > (U 2 ,l^) in Qr^ = T x (0,i22) and that there is 
a contact point either between Ui,U 2 at some xq G T, or between Vi,l^ at some 
(xo,yo) S ^r 2 (or both simultaneously). Proceeding exactly as in the proof of 
Proposition 13.31 we conclude that Ai < A 2 . Furthermore, since Vi ^ 1^ (thanks to 
the Dirichlet condition of V 2 at y = R 2 ), it even follows from the strong maximum 
principle that Ai < A 2 . 

Let us now prove the lower bounds on —A/j(a) (the upper bound was proved in 
Proposition 13.3p . Let Ua{x) > 0 be the principal eigenfunction associated with Aq. 
Letting z = one can check that Li{Ua^R,Va^R) = AR{a)Ua,R in (13.41) can be 
rewritten as 

-Dz" - 2D ia + ^\ z' + {-Da^ + K- AR{a))z = > 0. 

\ Ua I Ua 

We proceed by contradiction and assume that —Da^ + Aq — AR{a) < 0. Then any 
positive constant is a subsolution of the above equation satisfied by z. In particular, 
since z > 0 is periodic, we infer that it is identically equal to its minimum and thus 
constant. It then follows that 

n = {—Da^ + Aq — AR{a))z < 0, 

Ua 

which contradicts the positivity of Va,R. 
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It only remains to prove the other lower bound, namely 

-Aji(a) > da^ + /'(O) - d^. 

For fixed R > 0, denote uir = Arguing again by contradiction, if our lower bound 
does not hold then 


uj ;= 


da"^ + /'(O) + Aij(a) 


d 


>OJR>0. 


Averaging in x the equation for Va,R we see that ^R{y) := fj Va^R(x, y)dx solves 

^'kiv) + ^‘^^R{y) = 0 . 

Since ^R{y) > 0 for y G [0,R) and ^r{R) = 0 (recall that the eigenfunction 
Va,R{x,y) is positive up to the road and satisfies zero Dirichlet boundary condi¬ 
tions at y = ii), we get that 


^R{y) = C sm{uj{R - y)) 

for some constant C > 0. If w > lor then [0, R) contains at least half a period of 
$/£, hence ^r must have a non-positive value in [0, ii). This contradicts the strict 
positivity up to the road and the proposition is proved. □ 

We can now infer from the previous proposition that for any a fixed A/j(a) 
converges as R ^ -|-oo, and we thus define the generalized principal eigenvalue 

A(a) := lim Aij(a). (3-7) 

R —^-|-cx) 

It remains to prove that the associated pair of eigenfunctions also converges to a 
non-trivial limit as i? —>■ -|-oo. From now on, we will denote by (Ur, Vr) and Ar the 
principal eigenfunction pair and eigenvalue of the truncated problem (omitting the 
a dependence for convenience). 

Lemma 3.5. Assume that ||t7i?||L°°(T) = 1? then there exists a positive constant Ci, 
which is independent of R large enough, such that 

ll^it(')0)llL°°(T) > Cl- 

Proof. We argue by contradiction. If the proposition does not hold, then we can find 
a sequence —>■ oo such that 


Moreover, we assumed that ||C7fl||Loo(T) = 1- particular, for all k, there exists Xk 
such that 

\\UrJ\l°°{t) = U{xk) = 1 . 

Then the sequences of functions {URffjk and {VRf.{-,0))k are uniformly bounded in 
L°°(T). By standard elliptic estimates |17) . one can extract a subsequence such that 
UR^, —)• Uoo and x^ Xoo as k ^ -|-oo, where Uoo > 0 satisfies 

- 2DaU'^ + i-Da^ + y{x))Uoo = A(a)t/oo, 
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together with 


U 30 (Xoo) 1- 

By the strong maximum principle C/qo is positive, and by construction it is also 
periodic. By uniqueness of the principal eigenvalue Aq, of the periodic operator 

-DU" -2DaU' + ix(x)U, 


we get that 


A(o;) — — DcP' “t" Aq,. 


Because Aj:j(a) is decreasing with respect to i?, it follows that for all i? > 0: 


> —Da^ + Aq, 


which contradicts the previous proposition. □ 

Lemma 3.6. Normalizing with ||Vr(-, 0 )|| 2 ,oo|'’j') = 1, we have the upper estimate 

0 <VR{x,y) (x,y)eOij. 

Proof. Note that, if there exists some Rq such that Aj:jQ(a) + da‘^ + f'{0) < 0, then 
by monotonicity A/j(a) + da^ + f'{0) < 0 for all R > Rq. In this case, one can 
apply the maximum principle and conclude that Vr cannot reach a positive interior 
maximum. Therefore, recalling that 0 < VR\y=Q < 1 and VR\y=R = 0 < 1, we get 
that VR{x,y) < 1 and our statement holds. 

However, we only know that Aij(a) + da^ + f'{0) — d^ < 0 by Proposition 13.41 
Thus it remains to consider the case 

2 _ ^R{a)+da^ + /'(O) ^ ^ 

UJ — -;- > u 

d 

with 0 < uj < ojR := We denote again the average = Jr^VR{x,y)dx, which 

satisfies as before <I>^ + uj^^r = 0 and thus 

^it( 2 /) = Csm{uj{R - y)) 

for some C > 0. The fact that 0 < a; < ujr = tx/R of course ensures that ^R{y) is 
positive for y G [0, i?), but we also require that 

because of ddyV{x,0) = v{x)V{x,t)) — fi{x)U{x) < iyiV(x,0). Thus, 

—cj cotan (wii) < -J-. 

It is then easy to see that uj G [0,tJ], where uJ = uJ{R) is the unique solution of 

—tJcotan(a;i?) = -J- > 0 

in ( 0 ,^). 
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Introducing the function 


Vnix, y) := . sin {uj{R - y )), 

sin (wit j 

we claim now that this gives an upper barrier, that is 

V(x,y) G T X [0,i?] : Vr{x, y) < V R{y) . (3.8) 

Note first that Vr, Vr > 0 for y G [0, i?) and both are up to the boundary 
{y = R} where dyV r < 0, thus 

9* := min{0 > 0 : 9Vr > Vr for (x,y) G 0,r} > 1 

is well-defined. The inequality 9* > 1 immediately follows from our normalization 
on the road. 

We establish now (13.8|) . Assuming by contradiction that 9* > 1 and letting the 
elliptic operator L = —A — 2adx, then ^[1/^] = lo'^Vr and = oj‘^V r. As a 

consequence z := 9*Vr — Vr > 0 satisfies 

L[z] = co‘^9*Vr - u^Vr = (a;2 _ u:^)9*Vr + u}^{9*Vr - Vr) > 0 in Qr, 

because u <uJ and 2 : = 9*Vr — Vr > 0. 

Since Vr > Vr on the road and as 0* > 1 clearly z\y=Q > 0, so the strong 
maximum principle shows that 2 : > 0 for y G [0,i?). By the Hopf lemma and x- 
periodicity we also get dyz\y=R < 0. This easily implies that 2 ; > eV r in IIr for 
some small e > 0 and {9* — e)V r > Vr, which in turn contradicts the minimality of 
9* and thus entails our claim (13.81) . 

Recalling that 0 < w < lor = tt/R, we see that VR{y) = Csin(uJ{R — y)) is 
concave in y G [0, R\. Thus by definition of U: 

VR{x,y) < VR{y) < Vr(0) -b V'^(0)y = 1 - aJcotan(wii)y = l + ^y 

for all y G [0, i?] and x G T, and the proof is complete. □ 

Lemma 3.7. Normalizing with ||C^r||roo(t) = 1, we have 

ll^ij(-)0)llL°°(T) < C 2 
for some C 2 > 0 independent of R. 

Proof. Assume by contradiction that our statement does not hold. Suitably nor¬ 
malizing, we can therefore assume that ||Vr^. (•,0 )||ioo(']i’) = 1 and ||Lrj —)• 0 

for some subsequence Rj -|-oo. By (|3.7I) . Lemma l3.6[ and using standard elliptic 
estimates as before, we get up to extraction of a subsequence that Vr^ converges 
(locally uniformly on compact sets of fl) to a non-trivial function Voo > 0, which 
satisfies 


-dAVoo - 2dadxVoo - [da^ + /'(O) -b A(a)) Voo =0, x G T, y > 0, 
c^'9yVoo(x,0) = r(x)Voo(x,0), X G T. 


(3.9) 
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Moreover, standard elliptic estimates also allow to pass to the limit on the road 
equation and get 

i'{x)Voo{x,^) = 0 . 

Choose now any xq such that i'{xo) > 0. Then Kx)(xo,0) = 0, and by Hopf lemma 
dyVoo{xo,0) > 0. This contradicts the boundary condition in (13.9p and the proof is 
complete. □ 

Proof of Theorem \3.1[ We recall that A(a) was defined as the limit of Aj:j(a) as 
R —>■ +00 in (j3.7p . and postpone the proof of its concavity to the next subsection. 

Combining the above Lemmas l3.5ll3.6l and l3.7l together with elliptic estimates |17) . 
we can now pass to the limit as i? —)• +oo. We get a pair of non-negative and periodic 
eigenfunctions ([/„, V^) of (|3.3p in the whole field with A = A(a), which satisfy 

Ua{x)<l, Va{x,y) <C{l + y), 

after a suitable renormalization. The fact that these eigenfunctions are positive is 
a straightforward application of the strong maximum principle and Hopf lemma, as 
we already used extensively, and part (i) in Theorem 13.11 is proved. 

Let us now briefly check part (ii). Let some A be such that there exists a 
positive and periodic eigenfunction (U,V) of (j3.3p . Proceeding as in the proof of 
Proposition 13.41 (more precisely, the proof of the monotonicity of A/j with respect 
to R), one can check that A < AK(a) for any i? > 0. Passing to the limit R —)• -|-oo, 
it immediately follows that A < A(a). □ 


3.3 Further properties 

In order to complete the proof of Theorem 13.11 it only remains to prove the concavity 
of A(q:), which will play an important role in the study of the spreading speeds for 
exponentially decaying initial data. We exploit again the construction of A (a) via 
the more convenient framework of truncated problems. 

Proposition 3.8. The functions a i— >• A(a) and a e-)• AR(a), for all R > 0, are 
concave. 


Proof. Since the pointwise limit of concave functions is again a concave function, we 
only need to prove the concavity of A/j(a). Following the steps of |5l Proposition 5.7], 
we begin by showing that 


= sup min ^ ini-^-, mi 


(3.10) 


T 4> 

where 

T = {(<(),V') ^^^(T) X (C'2(L!fi)nC'i(n^)) : 

(^ > 0, V’ > 0 in T X [0, R), 'tp{y = i?) = 0 > dy'ip{y = R), E{4>, fj) > 0} . 

Note that the pair {Uo,,r,Vo,,r) G £ and, therefore, 

A ^ ^ / ■ i- c , L2i(f,'ip) 

Aij(a) < sup mm < ml-^-,mt 




T 4> if 
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Now proceed by contradiction and assume that there exists {(j), ip) € £ such that 


A_R(a) < min 


inf 


Li( 0 ,' 0 ) 


inf 


L2{(p, V’) 1 

V’ J ' 


Proceeding as before and thanks to the definition of the admissible set £, one can 
find some critical 0 > 0 such that Ua^R < 6(j) and < 9^p with either Uo,,r{xo) = 
^^(xo), Va,R{xo,yo) = 6'ip{xo,yo) or dyVa,R{xo,R) = 6dyip{xo,R) for some xq G T 
and 2/0 G [0, R). 

Note that {9(f) — U^^r, 9^p — Vol,r) > 0 satisfies 


' Ll{9(t) - Ua,R,9^P - Va,R) > Afi(a)(00 - Ua,R), 

L2{9ct> - U^^rM - V^,r) > hR{a){9'iP - V^^r), 
E{9(P-Ua,R,9p;-Va,R)>0, 

. {9P;-Va,R){y = R) = 0. 

By the strong maximum principle and the Hopf lemma as before, it follows that 
9(f) = Ua^R and 9)^ = which is a contradiction. Hence, (|3.10l) is proved. 

Let us now proceed to the proof of concavity. We first introduce 

= {(<^) : 3((/i>, Ip) & £ such that 0 = and p) = 


It is then clear that 

A / ^ fLi{(p,)p) . L2{(p,'ip) 

Aijfaj = sup mm < mt-^-,mt 


= sup min < inf 

{P>,' 4 >)€£c, \ ® 


(p ’ “djj p) 

—D^" + — up) 


inf 

lRx[0,R) 


—dip" — f'{0)pj 


Let ai > 0 and 02 > 0, and choose any {(pi,pJi) = e'^^^{(pi,pJi) G £ai, {4>2iP’2) = 
e"^^((/> 2 , ^^ 2 ) G £a 2 with {(pi,pii) G £. Fixing any t G (0,1) and defining a = tai + 
(1 — t)a 2 , we claim that 


Indeed since {(pi,pJi) G f we have that p)i,p )2 both vanish at y = R with non-zero 
slopes pi(x) ;= dyp)i{x,y = R), P 2 {x) := dyp) 2 {x,y = R), and it is then easy to 
check that p) = p)\p) 2 ~^ also vanishes at y = R with non-zero slope p = Pip\~^ and 
p) G The other conditions for {(p,p’) = {(p\(pl~^, p>\pjl~^) G £ also follow 

from straightforward computations, thus {(p,p)) = e“*(0, V’) G £a as claimed. 

Then one can check that 


. I. -Dp)" + yip - up) -dip" - f'{0)'^ 

mm < mt-^-, mf — 


Kx [0,iJ) 


p) 


^ ^ . \. f -Dip'l + pipi -uipi -dp)'l - /'(O)V’i ] 

> t mm < mt- - -^-, mf ' 


P>i 


, im -^- > 

Kx[0,i/) pJi J 

■d^2 - /'( 0)^2 1 


+ (l-«)mm ^- 

[ R p)2 IRx[0,R) Pj2 j 


As {(pi,pJi) and {(p 2 ,pJ 2 ) were chosen arbitrarily in respectively £a^, £ 02 , it follows 
that A/j(a) > tA/j(ai) -|- (1 — t)A/j(a 2 ). This concludes the proof. □ 


21 

















4 Spreading speed of solutions 


We are now in a position to prove Theorem 11.21 We begin by a characterization of 
the spreading speed c*{D) thanks to the generalized principal eigenvalue A(a) that 
we constructed in the previous section (see Theorem 13.11 above). 

Indeed by analogy with the single equation we introduce 

c*(T>) := minG [c|fPP,+oo), (4.1) 

a>0 a 

where cj^pp = 2\Jdf'{0) > 0. This c*{D) is well-defined thanks to the following 
inequalities 

max{Zla^ — Aq , da^ -|- f'{0)} < —A(a) < Ma, (4-2) 

which are in turn immediate consequences of Proposition 13.41 and (13.7p . Indeed note 
from (14.2p that —A has quadratic growth as a —)• -|-oo and —A(a) > f'{0) > 0, hence 
by continuity reaches its minimum as in (|4.1I) . Using again (|4.2I) . it is also 

clear that c*(D) > min = cVpp- 

In particular the equation 

A(q!) = —ca 

admits a solution a > 0 if and only if c > c*{D). Moreover, by the concavity of A(a) 
(Proposition 13.811 . if c > c*{D) then there are exactly two positive solutions. We can 
now establish the upper bound in our propagation result: 

Proof of the first part of Theorem \l.‘A For any solution (m, v) of (ll.3p - (ll.4p with non¬ 
negative and continuous compactly supported initial data we need to show that for 
any c > c*{D) and i? > 0 it holds 

lim sup {u{t,x) + v{t,x,y)) = 0. 

t^+oo x<—ct , 0<y<R 

Recalling the discussion at the beginning of Section O c > c*{D) is a necessary 
and sufficient condition for existence of a positive solution of the form (|3.2I1 to the 
linearized problem ()3.ip . Thanks to the KPP assumption f{u) < f'{0)u any such 
solution is also a supersolution of the original nonlinear problem (|1.3I) . 

More precisely, for any c > c*{D) choose some d G [c*{D),c) and a such that 
—A(a) = ad, and let (Ua,Va) be the associated positive generalized eigenfunctions 
from Theorem 13.11 The pair lies, up to multiplication by some 

constant, above the initial datum (uq,vq) at time t = 0. Recalling that the solution 
{u,v) is uniformly bounded, we can apply Proposition 12 .1 1 to get 

sup {u{t,x) + v{t,x,y)) < (||U„||lcx=(t)||U a||Loc(f^^)) . 

x<—ct , 0<y<R 

and the desired conclusion immediately follows. □ 

From now on, this section will be dedicated to the proof of the inner spreading 
estimate, namely the fact that for any 0 < c < c*{D) and R > 0, 

lim sup {\u{t,x) - U{x)\ + \v{t,x,y) - V{x,y)\) = 

t^+oo -ct<x<0 , 0<y<R 

where {U,V) is the unique positive and bounded stationary solution of (II.3p . 
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4.1 The construction of subsolutions 

In order to construct suitable subsolutions and obtain a lower estimate in the prop¬ 
agation Theorem 11.21 we first need the following technical result: 

Proposition 4.1. For fixed R> 0 the maps 

a !-)• Ania) £ C and a >->■ {Ua,R, Va^n) £ 

can be holomorphically extended to a complex neighborhood of the positive real axis 
{a > 0}. The complex-valued eigenfunctions still satisfy (|3.4I) together with some 
normalization, and moreover a Ua^R, Va,R are continuous with respect to the 

X M^) topology (here we mean the usual topology for the real 

and imaginary parts of each component U, V after identifying x 

Note that the restriction to truncated domains is important here. Indeed, the 
generalized principal eigenvalue A(a) in the infinite cylinder is not an analytic func¬ 
tion of a, see the comment after Proposition 15.11 later on. Since the equations are 
polynomial in a and the domain is bounded our statement follows by standard per¬ 
turbation theory combined with usual elliptic regularity and compactness arguments. 
We omit the details and refer to [181 Chapter 7], observing that the principal eigen¬ 
value Ajj(a) is isolated and has algebraic and geometric multiplicity 1, thanks to the 
uniqueness of the positive eigenfunction (see Proposition 13.31 above). 

For large i? > 0 we define 

c*R = min g (0, -^oo), 

a>0 OL 

which is the critical speed of the linearized problem in the truncated domain. Namely, 
system (13.ip (restricted to Qr) admits solutions of the exponential type (13.2p if and 
only if c > cfi. For the construction of subsolutions we shall need the second technical 
result below: 

Lemma 4.2. Let be the unique real and positive solution of — ^*R^*R- 

There exist some small r > 0 and 5 > 0 such that, for any c £ [c|j — S,cff), there 
exists a solution a{c) € C \ M of 

-A/j(a) = ac, 

which also satisfies \q*^ — Q:(c)| < r. 

Note that the existence and uniqueness of the positive solution of — Aj{(a) = 
follows from the definition of c|j, as well as from the concavity and analyticity of 
a I—>• —A/j(a). Then the conclusion of the lemma is an easy consequence of Rouche’s 
theorem and we omit the proof. 

We can now construct solutions of the linearized problem in some moving sets 
with speed c < c|j: 

Proposition 4.3. For all c £ [c|j — d,cfi), there exist real valued functions ui{t,x), 
vi{t,x,y) and h> 0 such that {ui,vi) is a solution of (13.ip in the moving set 

{{t,x,y) : t > 0, \xct\ < h and y £ {0,R)}, 
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and satisfies 


inf ui(t, —ct) > 0 , 
t>o^ ’ 


inf vi(t, —ct, R/2) > 0, 

t>0 ! / / 


and 


! ui{t,x)<0 ifx = —ct±h, 

vi{t, x,y) <0 if either x = —ct dszh or y = R. 

Proof. Let c G [c|j 5, c*jfi) and, by Lemma 14.21 there exists a solution 


a G c\: 


—A/j(a) = ac. 


(4.3) 

(4.4) 


By Proposition 14.11 there is a corresponding complex-valued eigenfunction pair 
{Ufi[x),VR{x,y)) of (13.4[) . By construction the functions 

{ui{t,x),vi{t,x,y)) = 3?e ([/^(x), Vr(x, y))^ 


are real valued solutions to (|3.ip on the truncated domain {{t,x,y) : 0 < y < R}, 
together with the Dirichlet boundary condition VR{-,y = R) = 0. 

Let us check that conditions (j4.3p and (14.4p are satisfied. Note first that, as 
c cjj, then a —)• a*^ and, by the continuity of {Ua,R, Va^R) with respect to a (see 
Proposition 14.11) , we can assume up to reducing 6 and without loss of generality that 

min3?e([/i{(x)) > ■imin[/^(x) > 0, min iR.e{VR{x,y)) > 0, (4.5) 

itET 2 xET xGT,0<y<R 

where {Ufi, Vfi) denotes the (normalized) principal eigenfunction pair of (13.41) with 

CX ^ ^R ■ 

Indeed for Ur the argument is straightforward, as well as for Vr away from the 
upper boundary y = R. More precisely, it is clear that for any 0 < e < R, then we 
can decrease 5 so that a is sufficiently close to and 

> I V^(x,y)> 0 . 

x£T,0<y<R—e A x£T,0<y<R—e 

In particular, ()4.3p is satisfied. 

In order to deal with y > R — e in (j4.5p . recall that the real valued function 
> 0 attains its minimum at any boundary point y = R, hence by Hopf lemma 
and periodicity, dyVfi^\y=R < -cq < 0. Using again Proposition 14.II and in particular 
the continuity with respect to a in the C'^’'’(n/j) topology (up to y = R), we see that 
there holds dy^e{VR)\y=R < —col2 < 0 up to reducing S again. Since VR\y=R = 0 
we get iR.e{VR) > 0 for y G [R — e,R), whence (|4.5p . 

It then immediately follows that, letting p = ?fiea > 0, u = Srma 7 ^ 0, and 
h = vr/|w|: 

{ui{t, —ct ± h),vi{t, —ct ± h, y)) = e^^^3?e (e^*^ {UR{—ct ± h), VR{—ct ± h, y))) 

< ( 0 , 0 ). 


Recalling that vi satisfies the Dirichlet boundary condition at y = R, this ends the 
proof of the proposition. □ 
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It is straightforward to check, by the uniqueness of the principal eigenvalue, that 
the function A/j(a) depends continuously on the parameter f'{0), in the L“^([ 0 , +oo)) 
topology. Thus, the speed also depends continuously on the parameter f'{0). In 
particular, for any small e > 0 , the proposition above still holds true if c|j is replaced 
by some which converges to as e —)■ 0, and f'{0) is replaced by /'(O)—e in (|3.1I) . 

Denoting by ui^s and the functions given by the above proposition applied 
to this perturbed problem, let then for any t > 0, \x + ct\ < h and 0 < y < R: 

u{t, x) = max(upe(f, x), 0 ) , u(t, x, y) = max(ui,£(t, x, y), 0 ). 

Then, by the KPP assumption and more specifically the regularity of / in a neigh¬ 
borhood of 0 , the pair {u,y) is, up to multiplication by some small constant, a 
generalized subsolution of the original nonlinear problem (II.3p in a compactly sup¬ 
ported set moving with speed c G — 5, c|jg). Moreover, choosing R large and e 
small enough, the speed c can clearly be chosen arbitrarily close to c*{D). 

More precisely, we have proved the following proposition: 

Proposition 4.4. There exist c G [0,c*) arbitrarily close to c*, and a pair {u,v) 
such that, for any 0 < k < 1, K,{u,y) is a continuous, non-negative and bounded 
(generalized) subsolution of (II.3p in a domain E x F where 

E := {{t, x) : —ct — h < x < —ct h} C (0, oo) x M, 


F := E X {0,R). 


Moreover, u{t,x -|- ct) and v{t,x -(- ct,y) are —periodic with respect to time, and 
satisfy that 

ini u{t,—ct) > 0 , ini y{t, —ct, R/2) > 0, (4-6) 


and 


uft, x) = 0 if X = —ct ± h, 

v{t, x,y) =0 if either x = —ct Eh or y = R. 


(4.7) 


Time periodicity and conditions (14.6p . (14.7p are immediate consequences of the 
construction of u, v above. While (14.611 ensures that the subsolution is not trivial, 
(HB is required in order to apply a comparison principle (after extending {u,v) by 
(0,0) outside of the strip |x -|- ct\ < h). 


4.2 Proof of the inner spreading theorem 

We can finally prove the inner spreading part of Theorem 11.21 Recall that (U,V) 
denotes the unique positive and bounded stationary solution of (|1.3I1 . Let us first 
prove the following two lemmas: 

Lemma 4.5. Let {u,v) he a solution of dLHD with bounded initial data (no,ro). 
Then 

limsup sup u{t,x) — U{x) < 0, 

t^oo xeM 

limsup sup v{t,x,y) — V{x,y) < 0. 

t—J-co a;GR,y>0 
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Proof. We begin by noting that, thanks to the KPP hypothesis, the pair 'y{U, V) is 
a supersolution of (II.3j) for any real number 7 > 1. In particular, the solution of 
the Cauchy problem (II.3p with initial data ( 7 C/, 7 I/) is non-increasing with respect 
to time. Moreover, by Proposition 12.11 it is bounded from below by the stationary 
solution (JJ,V). Therefore, it must converge as t —)• -|-oo to a bounded positive 
stationary solution, which must be {U,V) itself. 

On the other hand, as uq and vq are bounded and U, V have positive infimum, 
there exists some 7 > 1 such that {uq, vq) < 7(17, V). Applying again Proposition l2.ll 
we get the wanted conclusion. □ 

Lemma 4.6. Let [u, v) be a solution of don with bounded, non-negative and non¬ 
trivial initial data {uo,vq). Then {u,v) converges locally uniformly to {U,V) as 
t —)• + 00 . 

Proof. The previous lemma already proved that the supremum limits of u and v lie 
below U and V. Thus, we only have to prove that for any K > 0, then 

liminf inf u(t, x) — U(x) > 0, 

t^oo \x\<K,0<y<K 


liminf inf v(t,x,y) — V(x,y) > 0. 

t-i>oo \x\<K,0<y<K 

Recall that in a previous section, we have constructed arbitrarily small and compactly 
supported stationary subsolution pairs (U.{x),V_{x,y)) such that the associated solu¬ 
tion of the Cauchy problem (jl.3p converges locally uniformly to the unique bounded 
and positive stationary solution ([/, V) (see the proof of Lemma [2.4p . 

Moreover, as our system satisfies a strong parabolic maximum principle (see 
Proposition 12.II) . the solution {u,v) is positive everywhere for any positive time and 
we can assume without loss of generality that 


u{t = l,-)>U{-), n(t = l,-) >Z(-)- 


Applying again the comparison principle, the conclusion follows. □ 

We now go back to the proof of Theorem 11.21 Thanks to the two lemmas above, 
it only remains to prove that, for any 0 < c < c*{D) and i? > 0, then 

liminf inf u{t,x)— U{x) >t), 

t^oo -\-c*{D)>x>-ct 

liminf inf v{t, x,y) — V{x,y) > 0. 

t^oo —l—c*{D)>x>—ct,0<y<R 

From now on, c G (0, c*{D)) is fixed. Then, let c < c' < c*{D) close to c*{D) so that, 
by Proposition 14.4p there exists a moving subsolution {u,y) of (jl.3p with speed c'. 
Again, Proposition 12.11 implies that, for any t > 0, the solution {u,v) is positive in 
the whole domain. Thus, up to multiplication of {u,v) by some small enough k > 0, 
we can assume without loss of generality that for any x £ E and y £ F: 

u{l,x) > u{l,x) , v{l,x,y) > v{l,x,y). 

Furthermore, we also have that for any {t,x) £ dE (H ((l,oo) x M), 

u{t, x) > 0 = yft, x), 
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and for any {t, x, y) G dF n ((1, oo) x fi) 

v{t, x) > 0 = y_{t, X, y). 

Applying Proposition 12.21 (note that F D {y = 0} = E), one may immediately 
conclude that for any t > 1 and {x,y) G F: 

u{t,x) >u{t,x) , v{t,x) '>v{t,x). 

Now proceed by contradiction and assume that there exist 5 > 0 and some 
sequences ^ +oo, G [—1 — c*{D), —ctn] and G [0,i?], such that either 
u{tn,Xn) < U{xn) - S or v{tn,Xn,yn) < V{xn,yn) “ S. Note that, by parabolic 
estimates [2D], the sequences u{tn + s,a:„ + x) and v{tn + + x,y) converge 

locally uniformly to a solution {uoo,Voo) of which by Lemma 14.51 lies below 

{U,V). From the choice of and and a strong maximum principle argument, 
it is straightforward to check that (ttoo,^^oo) even lies strictly below {U,V). Thus, 
the sequence can be replaced by any sequence of points x'^ such that 
is bounded. In particular, we can assume without loss of generality that Xn = knL 
with kn € 7 j for any n G N. 

We can now let > 1. Then 

Xn + x) > Xn + x) = u(0, x), 

v{t'n,Xn + X,y)> Xn + x,y) = u(0, X, y), 

where the equalities follow from the time periodicity of the subsolution in the moving 
frame with speed c^ 

From Lemma 14.61 the solution of ([LSD with initial data iu{0,x),y{0,x,y)) con¬ 
verges locally uniformly to (JJ, V). Therefore, applying the comparison principle and 
noting that tn — ^ +oo as n —-|-oo, we get that {uoo,Voo) > {U, V) and reach a 

contradiction. This ends the proof of Theorem 11.21 


5 Speed enhancement by the road 

We have now computed the spreading speed c*{D) of solutions of (ll.3|l with com¬ 
pactly supported initial data, in the sense of Theorem 11.21 In this section, we will 
compare c*{D) with the spreading speed ^KPP of solutions of the single KPP equa¬ 
tion. In particular, we will prove Theorem 11.31 and show that the road accelerates 
the propagation (for compactly supported initial data) if and only if D > 2d. 

The proof will rely on Proposition 15.11 below, which compares A(q;) to the prin¬ 
cipal eigenvalue arising when looking for exponential solutions of the homogeneous 
KPP problem with no road. Furthermore, we will prove in the last subsection that 
A(a) also characterizes the spreading speed of solutions with exponentially decay¬ 
ing initial data. Therefore, Proposition 15.11 also infers whether such solutions are 
accelerated by the road or not. 


5.1 Eigenvalue enhancement by the road 


Noting that 


cTpp = min 
a>0 


da^ -|- f'{0) 
a 


da\pp + /(O) 

axpp 
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with axpp = Y "^-^7 the fact that c*{D) > d^pp if and only if D > 2d is a simple 
corollary of the following proposition: 

Proposition 5.1. 

(i) If D < d then 

Vo > 0 : —A(a) = da^ + f'{0). 

(ii) If D > d and a{D) = \J~^^ then 

VO < a < a{D) : —A(q:) = da^ + /'(O), 

Ma > a{D) : —A(q:) > da^ + f'{0). 

From {a) it is now clear that the generalized principal eigenvalue A(a) in gen¬ 
eral cannot be analytical in a, which contrasts with Proposition 14.11 in the case of 
truncated cylinders. 

Proof. Recalling from (|4.2I) that —A(a) > da^ + f'{0) we have for all a > 0 that 

- f'{0) - A(a) ^ ^ 
d “ ’ 

and that our statement amounts to determine whether w = 0 or w > 0. Letting [/ 
and V be the positive eigenfunctions of (13.3p from Theorem 13.11 integrating the field 
equation with respect to x, and denoting ^{y) = f^^jV(x,y)dx, it is easy to check 
as before that 

^"(y) - = 0. 

Since V(x, y) > 0 grows at most linearly in y (Theorem 13.ip this immediately implies 
that either d>(y) = ay + b for some a, 6 > 0 if w = 0, or <l>(y) = Ce~^y for some 
C > 0 if w > 0. Integrating now the boundary condition and the road equation, an 
explicit computation gives 

-d$'(0)= !{yU-uV{-,!))]dx = [Da^+ k{a)] [ Udx. (5.1) 

Jt Jt 

If either D < d, oi D > d and 0 < a < a{D), then Da^ < da^ -|- /'(O). From the 
above equality, we get 

-d$'(0) < [da^ + /'(O) + A(a)] [ U = -dw^ [ U. 

Jt Jt 

Thus <h'(0) > fjU > 0. It follows that oj = 0 (otherwise <I>(y) would be an 
exponential with slope —Cco < 0 at y = 0), which is exactly the desired conclusion 
in those two cases. 

In the last case D > d and a > a{D), then Da^ > da^ + f'{0) which, together 
with (15.ip . leads to 



As $'(0) > 0 when w = 0, it follows that a; > 0 and the proof is achieved. □ 



28 








5.2 Large diffusion limit 

Let us now end the proof of Theorem 11.31 by checking that 


0 < liminf 

D^OO y/D 


< limsup 

D^OO 


c*{D) 

Vd 


< + 00 . 


In order to obtain the upper bound we take first a = a{D) = y for D > d, and 
compute explicitly 

IW 

D-d' 

Because we chose exactly a = a(D) we are in case (ii) of Proposition 15.11 with 
—A(a) = da^ + f'{0), thus we just exhibited a solution a = a{D) of — A(a) = ac 
for this particular value of c. By definition of c*{D) this means that c*{D) < c, and 
it immediately follows that 

limsup^—< a//'(O) < + 00 . 

D^oo \D 

Turning now to the lower bound, recall from (I4.2jl that — A(q;) > max{Ila^ — 

Aq, + /'( 0 )}, where Aq S is the principal eigenvalue of+ 

/i(a:) on the torus. In particular we have 

—A(a) > p{a) := max{Ila^ — //i, + /^(O)}. 


da^ + /'(O) 
c :=- = U 


a 


Studying the piecewise polynomial and convex function p{a) for hxed H > d it is easy 

to check that min f y/D — d. Since c*(D) = min > min it 

«>o " “ V/'(o)+Mi ^ a>0 “ “ «>o “ 


immediately follows that 


lim inf 
D^oo 


c*{D) 

Vd 


> 


fjO) 

\/7^(o)'+7p 


> 0 . 


5.3 Exponentially decaying initial data 

Let us now turn to the proof of Proposition 11.41 By Proposition 15.11 we need to 
prove that c{a) = is the spreading speed of solutions for initial data with 

exponential decay of order e"^. This should be natural by now since for any hxed a 
and by construction of A(a) (see again Theorem l3.ip . c{a) is the smallest speed c such 
that there exists a positive solution of the linearized problem (|3.1I) of the type (|3.2I) . 

Proof of Proposition \l.f\ Let us hrst check the upper estimate. Note that although 
(17 q(x), Va{x,y)) is an obvious supersolution of (|1.3D . the difficulty lies in 
the fact that e^^Va{x, y) may decay as y —)• +oo and, therefore, it may not lie above 
the initial datum vq. 

Choose any c > c{a) and observe that ca > — A(a) by dehnition of c{a) = — 

An explicit computation shows that 

(I7i,Fi) := 
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satisfies 


dtUi — Dd^U 1 + fiU 1 — vVi\y=Q = (A(a) + ca)U i, t > 0,x G M, 

< dtVi - dAFi - /'(0)Fi = (A(a) + ca)Fi, t > 0, (x, y) G 0, 

—ddyVi\y=o = y,Ui — vVi\y=o, t>0 ,xGM. 

Also defining 

(1/2,F 2 ) := 

and recalling that —A(a) > da^ + f'{0), it is easy to check that 

I dtV 2 - dAV 2 - f'(0)V2 > (A(a) + ca)V 2 , t > 0, (x, y) G O, 

I -ddyV2\y=0 > fJ-U2 - l^V2\y=0, t > 0, X G M. 

Then, letting M large enough so that 

AI{A{a) + ca) min Ua > 1 ^ 1 , 

it is straightforward to check that M{Ui,Vi) + {U 2 ,V 2 ) is a supersolution of the 
linearized problem (13.Ij) and is moving with speed c > c{a). Up to multiplication 
by some large positive constant and from our assumptions in Proposition 11.41 it can 
be assumed to he above (tto,Xo) at time t = 0. It is then straightforward, applying 
Proposition 12.11 and noting that c could be chosen arbitrarily close to c(a), to reach 
the desired conclusion. 

Let us now focus on the inner spreading estimate. Proceeding as in the proof 
of Theorem [Q (see Section [4.21) . it is enough to find a non-trivial subsolution lying 
below (uq,vq) at time t = 0 and moving with speed c smaller but arbitrarily close 
to c{a). Let any R> 0 and define 

cr{q) = -^ g(Q,^ ag +00. 

a 

In particular Cfi{a) > c|j provided that R is large enough, and, by concavity of 
Aij(a), one has 

CR{a){a + ??) > -A/j(a + rf) 

for ry > 0 small enough. For simplicity we write below c = cr^o), and define 

if 0 < y < i?, 

where A and B are positive constants to be chosen later. 

One can check that for all f > 0 and x G M there holds 

dtU^ - DdlU^ + /iUi - vV^\y=o 

= A-Q- B [AflCa + ry) + (a + ry)c] < 0, 

and 

ddyV_l\y=Q T ^l^l|y=0 0 Ul|j/=iJ' 
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By regularity of / there exist a small sq > 0 and some 7 > 0 such that 

Vu G (—00, So] : f{v) > f'{0)v — max{0, 

Note that up to now f{v) was only defined for u > 0, hence we may indeed assume 
without loss of generality that f{v) = f'{0)v for u < 0 and, in particular, the above 
inequality holds. 

Choose first A > 0 small so that AVa^R < sq < 1, and then B > 0 large enough 
so that lJ_i < 0 and Y_i < 0 for all x + ct > 0. If moreover 7/ > 0 is small enough 
such that a + T] < a(l + r) then for x < —ct it is easy to check that 

max{0,7Zi}'+" < 


whence 

dtVi - dAZi - /(Zi) = [/'(0)Zi - /(Zi)] + [dtVi - dAZi - /'(0)Zi 

< max{0,7Zi}^+’' + U • 0 - + r?) + (a + 


< Q[ot+-n){,x+ct) 


{lAVa^Rf^"' - B{AR{a + ?/) + (« + T])c)Va+r,,R 


< 0 


if i? > 0 is large enough. On the other hand when x > —ct the same inequality 
dty_i — dAV_i — f{y_i) < 0 easily follows from the fact that y_i < 0 and thus 
/(Zi) = /^(0)Zi- Therefore the pair 


U_ := maxlO,^^}, 


f max{0,Zi} if 0 < y < ii, 

I 0 elsewhere. 


is a generalized subsolution of (II. 3p in the sense of Proposition 12.21 By construc¬ 
tion this subsolution moves with speed c = c/j(a) = — —which can be made 
arbitrarily close to c(a) = —by choosing R large enough. Moreover, with our 
assumption on the initial data and up to multiplication by a small constant, this sub¬ 
solution can be assumed to lie below (uq, uq) at time t = 0 (observe that U_ = y_ = 0 
for X > —ct). As mentioned above, one may then proceed as in Section [4.21 to end 
the proof of Proposition 11.41 □ 


6 Some extensions 

The crucial feature of the road-field system is that it is of a cooperative type, which 
in particular allowed us to construct and characterize the spreading speed through 
a family of principal eigenvalues and Krein-Rutman theory. We briefly discuss here 
some extensions to more general frameworks where a similar approach could be 
performed. 

Truncated fields: We first highlight here the fact that the generalized principal 
eigenvalues A(a) were obtained as the limits of the principal eigenvalues Aij(a) of 
some truncated field problems (see Section [3] for the details). Therefore, it should be 
clear from our proofs that a similar result as Theorem 11.21 holds for the restriction 
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of (11.31) to a strip domain {{x,y) G M? : 0 < y < i?} with a Dirichlet boundary 
condition v{y = R) = 0. 

More precisely, provided that R is large enough, the solutions spread with the 
speed G (0, oo) along the road, where was defined in Section |4l When y and v 
are positive constants, this was proved in |26| . along with the fact that there exists 
a critical diffusion D above which c|j > c^pp and below which Cp < Cppp. This 
relied on the monotonicity of Cp with respect to D, which remains an open problem 
in our periodic framework. 

Noting that the propagation is slowed down by the Dirichlet condition on the 
upper boundary and in order to isolate the effect of the line of fast diffusion, one 
may also want to compare Cp with 


Ckpp,r ■— 2 



d^TT^ da^ + /'(O) - df^ 

—= mm-^ 

AR^ Q>o a 


which one may check to be the critical spreading speed of solutions of 


dtv - dAv = fiv), 


t > 0, X G M, y G (0, R), 


dyv{t, X, 0) = 0 = v{t, x,R), t > 0, X G 


( 6 . 1 ) 


Because the main idea of the proof of Proposition 15.11 (thus of Theorem 11.31 and 
Proposition 11.41) roughly consists in determining the direction of the (average) flow 
between the road and the field, the above problem naturally arises. Note also that, 
in the limit case R = +oo, it is clearly equivalent to dEB. 

Therefore, by a similar argument as in Section [5l one may prove that 


—A/j(a) < da^ + f'{0) — d 


vr 




if 


a < aR{D) := 


'no)-d 


4iP 


D-d 


while the opposite strict inequality holds if a > aji{D). 

Then, letting ap be the unique (thanks to the analyticity of —Ar) positive 
solution of CpO = — A/j(a), it is straightforward to obtain the analogous of Propo¬ 
sition 11.41 in the truncated field case. Namely, solutions with exponential decay of 
order e“*, where a < ap, are accelerated by the road in the sense that they spread 
with some speed 

-AR{a) ^ da^ + /'(O) - d-^ 


CR{a) = 


> 


t2 

''47F 


a a 

if and only if a > aR{D). The right-hand term of the above inequality is, of course, 
the spreading speed of solutions of (16.ip with decay of the same order e“^. 

However, it is not clear in general whether the minimal speed Cp is larger than 
^*KPP,R- This is related to the fact that we cannot always locate ap where the 
minimum in the definition of Cp is reached. In particular, the above discussion 
does not provide a precise diffusion threshold for the acceleration of solutions with 
compactly supported initial data (or more generally, initial data that decay faster 
It is only straightforward to infer from the above that Cp < Cpppp 


ii D < 2d as well as, using Proposition 13.41 and more specifically the estimates on 
—AR{a), that Cp > c*j^ppR for large enough D. 
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Including KPP reaction on the road: The key assumption, which in particular 
we used extensively in the investigation of the effect of the road, is the fact that the 
heterogeneity is limited to the exchange terms. Therefore, one may easily check that 
our method still applies when replacing the road equation in fll.3jl by 

dtu — Dd^u = u{x)v{t, X, 0) — ^{x)u + g{u), 


where g{u) satisfies a KPP type assumption 


g{0) =0, 


aju) 

u 


is non-increasing, 


and 


9{u) 

u 


< 0 for large enough u. 


We leave the details of the construction of the spreading speed to the readers: the 
proof of Section [3] applies to the letter up to minor changes (by simply adding 
g'{0)U to the operator Li, g'{0) in the definition of M„, and denoting by Aq, the 
principal eigenvalue of -D-^ — 2aD-^-\-g, — g'[Q)). For the sake of simplicity, let us 
just assume that we have constructed for this new problem a family of eigenvalues 
—Ag(a) > da^ -\- f'{0) and a spreading speed 


1(D) = min 
' «>o 


-Agia) 


a 


We now focus on Proposition 15.11 The same straightforward computation leads to 
the critical equality Da^ + g'{0) = da^ + or equivalently a = ^ ^ 

which separates whether —Ag(a) is larger or equal to da‘^ + f'{0). Using the fact that 

^*KPP ~ rninoo and that this latter minimum is reached for a = \J ^ , 

one then immediately gets that c*(D) > c*^pp if and only if 


D > 2d 


d 


am 

f'ioy 


This extends a result of m in the periodic exchange framework. Proposition 11.41 
also naturally extends and the solutions with exponential decay of order e"^ at time 
t = 0 are accelerated by the road if and only if Da^ + g'{0) > da^ + 


The general periodic framework: In the above extensions and as we already 
pointed out, we restricted ourselves to problems where only the exchange terms 
depend on x. The main reason is that we can only give an accurate diffusion threshold 
for the acceleration phenomena in such a framework (see the averaging argument in 
Section [5]). On the other hand, the construction of principal eigenvalues in the 
truncated field mostly relies on the classical Krein-Rutman theory, thanks to the 
monotone feature of our system. 

This means that this part of our argument, along with the existence of a spreading 
speed (at least in the truncated problem) could possibly be extended to a much more 
general framework as long as the x-periodicity is maintained. This would include 
heterogeneous reaction or diffusion (in both x and y) and, perhaps more interestingly, 
curved roads. 

Nonetheless, to extend Theorem 11.31 would then be a completely open but inter¬ 
esting problem. It of course seems unlikely that there always exists a simple formula 
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for the diffusion threshold. To exhibit such a threshold, one may need to prove 
the monotonicity of the principal eigenvalue on some diffusion parameter, which we 
avoided here. More generally, we hope to be able in a future work to further charac¬ 
terize this eigenvalue and, in a similar fashion as for the single equation |22| . derive 
not only properties on the dependence of the spreading speed on diffusion, but also 
on various parameters and the exchange terms. 
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